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Abstract. A zero-gap state (ZGS) has been found in a bulk system of two-dimensional organic conductor, 
a-(BEDT-TTF)2l3 salt which consists of four sites of donor molecules in a unit cell. In the present paper, 
the characteristic of the ZGS is analyzed in detail and the electronic properties are examined in the 
vicinity of the Dirac point where the conduction and valence bands degenerate to form the zero-gap. The 
eigenvectors of the energy band have four components of respective sites, where two of them correspond to 
inequivalent sites and the other two correspond to equivalent sites. It is shown that the former exhibits an 
exotic momentum dependence around the contact point and the latter shows almost a constant dependence. 
The density of states of each site close to the Dirac point is calculated to demonstrate the temperature 
dependence of the local magnetic susceptibility and the local nuclear magnetic relaxation rate. Further, 
the robust property of the ZGS against the anion potential is also shown by using the second-order 
perturbation. 

PACS. 71.10.-w Theories and models of many-electron systems - 72.80.Le Polymers; organic compounds 
(including organic semiconductors) 



1 Introduction 

Organic conductor, BEDT-TTF 

(bis(ethylene)dithiotetrathiafulvalene) salt has been stud- 
ied extensively since the various kinds of the transfer ener- 
gies between BEDT-TTF molecules show exotic phenom- 
• ena in condensed matter physics pQ. Recently, a quasi- 
two-dimensional conductor, a-(BEDT-TTF) 2 I 3 salt pp]. 
becomes remarkable due to the theoretical finding of a 
massless Dirac particle under pressures 4 , 5 based on the 
experimental data of transfer energies [6] . Such a particle 
exhibits the the zero-gap state (ZGS) where the valence 
band and the conduction band touch each other at a mo- 
menta, ±fco in the Brillouin zone, and a Dirac point is 
produced. As the noticeable properties due to transfer en- 
ergies, the location of fco varies under pressure and the 
energy band around the Fermi energy is described by a 
linear dispersion but with an anisotropic Fermi velocity, 
i.e., by the tilted Weyl equation [7]. 

The intriguing property in the a-(BEDT-TTF) 2 I 3 salt 
appears in the transport phenomena, which come from a 
zero-gap conductor. The Hall coefficient decreases rapidly 
with a power law as a function of temperature [5], while 
the electronic conductivity, a, stays almost constant (a ~ 
h/e 2 per BEDT-TTF layer) [5]. Such a behavior of the 
transport property has been successfully explained in terms 
of the ZGS PEED. 



We note that there are several materials which show 
the Dirac fermions. The Dirac particle in condensed mat- 
ter physics was first discovered in graphene (single layer 
graphite) [i~2lll3) . The graphene shows the isotropic lin- 
ear dispersion around the corner of the Brillouin zone pU 
[T5] , The structure of the Landau level [TB1IT7] . anoma- 
lous transport phenomena [T8lfT9ll20[l21[[25] , and impurity 
effect p3 l l24 l [25 | [26 l [27 ] , of the graphene have been pro- 
posed. Further, the bismuth has a small band gap with 
anisotropic velocities, and is described by the Dirac equa- 
tion [28 . The large diamagnetism appears due to the inter- 
band effect of magnetic field [35] . The ZGS also occurs in 
the HgTe quantum well by changing the thickness of the 

weii [ig. 

There are still several issues on the ZGS of the a- 
(BEDT-TTF) 2 I 3 salt, which are not yet clarified, (i) It is 
not obvious how the ZGS is constructed by BEDT-TTF 
molecules with several transfer energies. There are four 
molecules located at A, A', B, and C sites in a unit cell, 
where A and A' are equivalent sites due to the inversion 
symmetry. The features of these three kinds of molecules 
(A, B, and C) are different with each other. There ex- 
ists the charge disproportionation even at high tempera- 
tures [3TII32"] . and the relation is found as (nc) > {tia) = 
(tia>) > (tib) [33], where (n a ) denotes the amount of the 
charge. The difference in the magnitude of (n a ) becomes 
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large under pressures [34] . Further, it has been found re- 
cently that the local magnetic susceptibility Xa shows a 
relation, \c > Xa — XA' > Xb [34], which has the same 
relation as that of the charge disproportionation. In the 
present study, we examine the respective contribution of 
four different molecules to the ZGS, and show that the 
ZGS does not imply similarity between the susceptibili- 
ties and the charge disproportionation. 

(ii) It is known that the ZGS is robust against pressure 
although the location of the contact point, ±ko, varies 
due to variation of transfer energies under the pressure. 
Recently, the ZGS is examined in the presence of the anion 
potential, which gives the local potential on conduction 
electrons at each molecules [35 . It is of interest to analyze 
the effect of the anion potential on the stability of the ZGS 
and the location of the contact point. 

In the present paper, we study the role of the respective 
sites on the ZGS by calculating the local electronic state 
in the vicinity of the Dirac point. In SJ21 the energy bands 
and charge disproportionation are calculated for the a- 
(BEDT-TTF)2l3 salt by using transfer energies estimated 
from the first principle calculation. [36 The perturbation 
theory is also applied for the analysis of the electronic 
state close to the Dirac point and anion potential. In £j3l 
the electronic state close to the Dirac point is calculated by 
focusing on the local properties of four donor molecules. In 
SJU the effect of the anion potential is examined. In §5\ the 
local density of states is calculated and is applied to the 
estimation of the temperature dependences of magnetic 
susceptibility and the nuclear magnetic relaxation rate. 
Finally, the summary is given in fJS] 



2 Band calculation 




**x(b) 

Fig. 1. Structure of the conducting plane of the q-(BEDT- 
TTF) 2 I 3 salt. The ellipse denotes a BEDT-TTF molecule and 
the unit cell of the salt is drawn by the gray square. The bonds, 
al, a2, a3, 61, 62, 63, and 64, denote the nearest neighbor elec- 
tron hoppings and al', a3', and a4' denote the next nearest 
neighbor ones. Anions exist above and below dashed rectan- 
gles of conducting plane. 



Hubbard model, 

H = tiu:jl3c\ a! jCjl3e + J]] Iac\ acr Ciaa 

(ia,j/3),a ia.tr 

ia. 

+ V ^iP C \ a a C )f}*i C iP°' C iw (!) 

n.7i. : aa' 



The conducting layer of a-(BEDT-TTF) 2 I 3 is shown in 
Figure [U which are stacked with anion layers alternately. 
The notations of the bonds, al, a2, a3, 61, 62, 63, and 64 
in Figure [TJ are the same as those in Ref. [3], and quanti- 
ties al', a3', and a4' denote those at next nearest neigh- 
bor sites, which were introduced by Kino and Miyazaki 
[36] . In the present calculation, we adopt transfer energies 
evaluated from the first principle calculation at T = 8 K 
under the ambient pressure [3B]. Those in the unit of eV 
are given as t al = -0.0267, t a2 = -0.0511, t a3 = 0.0323, 
t bl = 0.1241, t b2 = 0.1296, t b3 = 0.0513, t b4 = 0.0152, 
t aV = 0.0119, t a y = 0.0046, and t aV = 0.0060. The unit 
of the energy is taken as eV, hereafter. Other choice of 
transfer energies is discussed later. 

The charge disproportionation in the a-type BEDT- 
TTF salt (Fig. [I]) is examined by considering the extended 



where c iaa denotes the annihilation operator of the fcrmion 
at the a-th site in the i-th unit cell with the spin, <r(=f , J.). 
In Figure [TJ sites A, A', B, and C correspond to 1, 2, 
3, and 4, respectively, and cv,/3 = 1, 2, 3, and 4. The 
quantity ti a: jp is the transfer energy between nearest or 
next nearest sites. U is the on-site Coulomb interaction. 
Via-.jp is the off-site Coulomb interaction between (i, a) 
and (j,/3) sites, where the interaction is taken as V a for 
al, a2, and a3 bonds, and V b for 61, 62, 63, and 64 bonds. 
The local potentials at A(A'), B, and C sites are given 
by ii(= ^2), ^3, and I4, respectively, which are examined 
in fJU The potential comes from the anion, I3, where the 
valence is given by —1 due to an electron transfer from 
two BEDT-TTF molecules. The anions also form a layer 
which is located between conducting layers. The position 
projected into the conducting plane (layer) is shown by 
dashed rectangles in Figure 1 [2]. 
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We treat interactions within the mean-field approxi- 
mation given by 

C iaa C jficr' C jP<?' C iota 

-> (n a <r)& j/Sa ,c j0<x < + c\ aa c iaa (nfj ol ) - (n aa )(n 0a ,) , (2) 

where we take into account only Hartree terms and discard 
the exchange term. The quantity, {(y iaa Ci aa ) — (n aa ) is 
calculated self-consistently. Using the Fourier transform, 
CiQtr = (1/y/N) J2k c kaae lk ri , eq. dTJ> is rewritten as 



Hmf = ^(4.. <2„ 4 4 J ( T (k) + I) 

k.a 




(3) 



where T(fe) = [t a p(k)], and I — [I a S a p\. The 4x4 matrix 
elements, t a p{k) (a, (3 = 1, • • • ,4), are given by 



in(fc) = 2t i' cos k y + 
tu(k) = t b2 +t b3 e ik * , 
*22(fe) = 2i i/ cos ky + 
t2i(k) = tu + 4ie lfex , 

t a p{k) = tp a {k)* , 
with 



C 2 



tu(fc) 
«23(fc) 

*33(fc) = 2ia3' COS fc y 
£44 (fc) = 2i a 4' cos k y 



toe 1 * 



t b2 e i{k * +ky) 
Cs, 

C4 , 

(4) 



d = C/(m)/2 + 2V a (n 2 ) + 2V h ((n 3 ) + (n 4 )), 
C 2 = U(n 2 )/2 + 2V a (m) + 2V b {(n 3 ) + (n 4 )), 
C 3 = U(n 3 )/2 + 2V a {r H ) + 2V b ({ni) + (n 2 )), 
d = U(m)/2 + 2V a {n 3 ) + 2V b ((n 1 ) + (n 2 )), 



(5) 



where (n a ) = + (n a i). In eq. the wave num- 

bers k x and k y are measured by 1/a, where a denotes a 
lattice constant. Note that, in the ZGS, there is no mag- 
netic moment at any site i.e., (m a )(= (n a -\) — ( n ai)) = 0. 
The relation, (ni) = (n 2 ), is satisfied due to the inversion 
symmetry. Quantities (n a ) and the chemical potential /x 
are calculated self-consistently by 

£(T(fc) +1- »E)d 7 (k) = £ 7 (fc)d 7 (fc) , (6) 
fe 

K) = ^Ei^(fc)i 2 /(c 7 (fe)), (7) 



where f(x)(= l/(e^ _M ^ T + 1)) is the Fermi distribution 
function, and T is a temperature (The unit of T is given 
as eV, and Boltzmann constant is taken as unity). For 
simplicity, we set T — and /(£ 7 (fc)) in eq. (JSj> is replaced 
by a step function, 6(ef — £ 7 (fe)). E is the unit matrix. 
The a-th component of the vector d 7 (fc) (a = !,••• ,4) 




Fig. 2. Energy band, £ 7 (fe) (7 = 1,2), showing Dirac cone 
in the vicinity of the contact point, feo = (0.686-7T, — 0.4437r). 
The Coulomb interactions, U = 0.4, V a = 0.17, V b = 0.05, are 
taken into account in addition to the transfer energies, which 
are estimated by the first principle calculation at 8 K under 
ambient pressure. We set fi(fco) = £2(fco) = 0. 



is defined as d aj (k). The energy band, £ 7 (fe) (£i(fc) — 
&(k) > £,3(k) > £t(fe)), and the eigenvectors <f 7 (fc) are 
obtained by diagonalizing eq. ©. 

We choose parameters, U, V a , and V b , in order to ob- 
tain the observed (n a ) [33], in which the largest electron 
number is given at the C site and the smallest one is given 
at the B site. In the present calculation, these parame- 
ters are taken as U = 0.4, V a = 0.17, V b — 0.05, and 
I a = 0.0 [33], which give the ZGS with contact points, 
fc = ±(0.686vr, -0.443tt). Notice that the case U = V a = 
V b — also gives the ZGS as shown in Ref. [36]. For the 
present case, the energy band of the Dirac cone around 
the contact point is shown in Figure [2] The charge dispro- 
portionation is given by (m) = (n 2 ) = 1.460 for A and A' 
sites, (n 3 ) — 1.369 for B site, and (714) = 1.711 for C site. 

Next, we examine the electronic state close to fco within 
the perturbative method. Instead of the band index 7, we 
use £ = 1, 2 for valence and conduction bands and 77 = 3, 4 
for the lower two bands. The Coulomb interaction and the 
charge disproportionation are determined unperturbative 
method, i.e., within eqs. ©-([H]). We construct the effective 
Hamiltonian as follows: 



H eff =H Q + V + V 



1 



En — Hn 



-V ■ 



(9) 



where V is the perturbative Hamiltonian and the explicit 
form is given in the next section. Hq denotes the Hamil- 
tonian, eq. ([3]), at fco, with the effect of the charge dis- 
proportionation ((rii) = (n 2 ) — 1.460, (n 3 ) — 1.369, 
(at) = 1.711, for U = 0.4, V a = 0.17, V b = 0.05, and 
1 = 0). Since the system exhibits the ZGS, the upper 
two eigenvalues of Hq degenerate on the Fermi energy, i.e. 
Eq = 61(^0) = £2 (fco) = 0, and the lower ones are given 
as £ 3 = £ 3 (fc ) = -0.218, and E 4 = £ A (k ) = -0.404. 
The eigenstates of H described by the Luttinger-Kohn 
representation [37] for £^(fco) and ^(fco) are given by |£) 



4 



S. Katayama el. at: Electronic Properties Close to Dirac Cone in «-(BEDT-TTF)2i3 



(C = 1,2) and \rf) (77 = 3,4), respectively. Note that \Q 
can not be determined uniquely due to the degeneracy. 
The detail forms of d\ and d 2 are discussed in the next 
section. On the other hand, |r/) can be described uniquely. 
Actually, the eigenvectors, d 3 (k ) and d±(k ) are calcu- 
lated as 



d 4 (k Q ) = 



0.5154 \ 
-0.1944- 0.4477i | 
-0.0944 + 0.4392i 
-0.2410- 0.4571i/ 

0.5190 \ 
0.3211 + 0.4077i | 
-0.4332- 0.031H 
-0.4233 + 0.3057i/ 



(10) 



When one of the linear combination of \ip^) — ai^(k)\l)- 
a2^(fc)|2) is considered, we obtain the equation for the en- 
ergy (AE^(k)) up to the second order of V given as, 



E 

£"=1,2 



(civic") + E 
77=3,4 



(C\v\v)(v\v\C") 

6? Oo) 



a C"C'( fc ) 



= AE c (k) a<:c (k). 



(11) 



In the next section, we discuss the ZGS using both the 
exact and perturbative calculations. 



(a) 



(b) 




Fig. 3. ^-dependence of the absolute values of the eigenvec- 
tors, |d«x(0)| (a) and \d a2 (6)\ (b) (a = A,A',B,C). The open 
circles are obtained from eq. (I21|) . The radius S'(= |fc — fco|) 
is set as 5' = 10 _7 7r. 8 denotes the angle between fe — feo and 
kx~ cixis. 



equation is written as 
fe 



«21 *22 



a 2C (fc) 



=fc 



W + «3 «X - 1^2 
«X + ™ 2 U - « 3 



a 1( (k) 
a 2C (fc) 



: E fe ' 

p=0 



a iC (fe) 
* a 2C (fc) 



a 2C (fc) 



(13) 



3 Electronic state in the vicinity of the Dirac 
point 



The velocity, i> 77 ' (7, 7' = 1, • • • , 4) is given by 



*77' = E^*7 ^f( fe o)^/37'^ 

ap 



(14) 



In g3] and SI the Dirac point fc = (0.686tt, -0.443tt) is 
taken as the origin of the wave vector (fe — feo — > fe). 

First, we discuss the eigenvector with k close to feo- 
FigureOshows the ^-dependences of |d Q i(0)| (a) and |<i a2 (#)| 
(b), which are calculated by diagonalizing eq. ijHJ) ((k x) k y ) = 
6' (cos 6, sin#)). The notable feature of Figure [3] (a) is the 
existence of a node for the components of both B site 
and C site, e.g., ^31(0) = and dn(7r) = 0. For ~ 0, 
|d4i (0)| is the largest, whereas |d 3 i(#)| becomes the largest 
for 9 ~ 7r. On the other hand, |d a i(#)| for a = 1,2 has 
weak ^-dependence. It is also seen the relation, |d a i(0)| — 
\d a2 (9 + 

The behavior of Figure[3]is analyzed by the tilted Weyl 
equation. For small |fe|, the Hamiltonian, V, is written as 



V(k) = \ko,a) 



dT 
fe-^(feo) 



(12) 



■4) 



where |fe,a) is the electronic state at site a(= 1. 
with the wave number fe. By substituting eq. (fT2"|) into eq 
(fTT|) and discarding the second order term, the tilted Weyl 



where d ai is defined as d a Q — d a ^(k' ) (£ = 1,2) and 
(feo) (77 = 3,4). Since the components of the 
eigenvectors of the conduction (d a \) and valence (d a2 ) 
bands can not be determined at feo, the vector k' = 
(6" cos 9, 8" sin 9) is introduced, where 5" is an infinitesi- 
mally small quantity and 9 is the angle between k' and 
fc x -axis. Velocities, vq, tlx, v 2 , and V3, are expressed as 
Vq = (vu + 6 22 )/2, v 1 = Re(v 12 ), v 2 = -lm(v 12 ), and 
v 3 = (^11 — v 22 )/2, and 00 and cr; (i = 1, 2, 3) denote unit 
matrix and Pauli matrix, respectively. The eigenvalue of 
eq. (fT3|) is calculated as 



AEi, 2 (fc) = k-v ± 



\ 



P =i 



--v 0x k x + v 0y ky ± ^jv^kl 



v-ki 



(15) 



and 



With V% = + + 7j| x , = + ^ t^J,, 

V xy = v\ x v\y + v 2x v 2y + V3 X V3 y . It should be noted that 
i>i and v 2 depend on the wave number fc , whereas the 
velocity vectors v and v 3 and quantities, V xl Vy, and 
V X y are independent of the choice of the base. 
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We choose d aj as the eigenvector of t a p(ko) with 
0, i.e., 



Xi = di 



x 2 = d 2 = 



0.5194 
-0.4429 + 0.2713i 


0.6623 + 0.1473i 

0.4419 
0.3703- 0.241H 
0.6125- 0.484H 



(16) 



The velocity vectors are given as 



v a = (-0.0389, 0.0048), v 1 = (0, -0.0005), 

v 2 = (0, 0.0439), v 3 = (0.0515, 0.0009). (17) 

These vectors, are drawn in Figure |H In this case, the 
effective Hamiltonian can be approximately written as 



i _ fvo x k x + v 0y ky + vk x 



H' = 



— lv'ky 

voxkx + voyky - vk x 



(18) 



where v 0x = -0.0389, v Qy = 0.0048, v' = v 2y = 0.0439, 
and v — vsx = 0.0515. Note that the effective Hamiltonian 
of eq. (fT8| differs from that of Ref. [10], since the choice 
of k' Q , i.e. the bases at fco, are not the same. However, the 
result does not depend on such a choice. The eigenvalue 
{AEq) and eigenvector (a^<) of eq. fig]) are given as 



AEi = (j = v 0x k x + v 0y ky + \jv 2 k x + v' 2 k 2 , 

AE 2 =J = v Qx k x + voyky - yjv 2 k% + v' 2 k 2 , (19) 

a n {0)\ _ /cos(0/2) \ 
a 21 {6) J \iaa(0/2)J' 

a 12 (8) \ _ ( sin(0/2) \ 
a 22 (6) I -i oos(e/2) J ' 



(20) 



where the relation v ~ v' is used for the calculation of 
d££/(fc), for simplicity. From eqs. (Til)]) and (|2H)l . the 6*- 
dependence of d aj (7 = 1,2) is obtained as 



di(0) 



a 2 i(0) a 22 (0) 



(21) 



o 



0.06 




uj=0. 


0.04 




.-*"' '' ■'/( 1 


0.02 


&'' 


v 2 











u (I c ) 


i ©-*-•-- - 



-0.04 



0.04 

x-component 



Fig. 4. Velocity vectors, v p (eq. (|17[) ) and u p (eq. (|30|l ). where 
the bases are given by eq. (|16[1 . b p are shown by arrows. For 
I a = -Tb = 0, and —0.2 < Ic < 0.2, w p (p = 0, 2, 3) are shown 
by dotted line where open circles correspond to Ic = 0.2. u\ 
is small compared with other velocities (|«i| ~ 10~ 3 ). 



of an ellipse as shown in Figure [5j The angle, </>, between 
the long axis of the ellipse and fc^-axis, and the ratio, a/6, 
between the long and short axes are calculated as 



= — tan 
2 



C3 



ci - c 2 



a 
b 



' ci + c 2 + ^/(ci - c 2 ) 2 
ci + c 2 - ^/(ci -c 2 ) 2 



(23) 



respectively. In the present case, one obtains <j) < and 
a > b. Further, by introducing a polar coordinate, k x — 
rcos9 and k y = rsin^, eq. (|22|) is rewritten as 



V4A 



1 f ^ - A 2 



2A, 



= F(0)w, 



(24) 



This gives |d 4 i (0)| = ^32(0) I - 0.678 1 cos(0/2)| and |d 3 i(0)| 
|c?4 2 (0)| ~ 0.785| sin(0/2)| which reproduce well the exact 
one as shown in Figure [3] (open circle) within the accuracy 
of -0.03. 

Next, we discuss the detail structure of the Dirac cone. 
Equation (fT9|) is rewritten as 



c\k x + c 2 k y + c 3 k x k y + c A k x + c 5 k y = 1, (22) 

where the coefficients are given by c\ = (v 2 — Vq x )/uj 2 , 
c 2 = (v' -vl v )/uj 2 , c 3 = 2v 0x v 0y /uj 2 , c 4 = -2v 0x /uj and 
C5 = 2vo y /uJ. Equation (|22p represents the general form 



where A\ = ci cos 2 + c 2 sin + C3 cos sin and A 2 = 
c 4 cos0 -I- c 5 sin 0. For U = 0.4, V a = 0.17, V b = 0.05, the 
0-dependence of F(6) is drawn in Figure[5l where parame- 
ters are given as a/oj = 48.00, b/co = 35.20, X /ui = 36.09, 
Yq/uj = -7.082, and = -0.089tt ((X ,Y ) denotes the 
center coordinate of the ellipse.). The general form in 
terms of Xo and Yq is shown in Appendix. 

Let us compare the present result with that of the 
graphene. For the latter case, the Hamiltonian of eq. (Ti"3")) 
is rewritten as H' = v(k x <j\ + k y a 2 ), which denotes the 
isotropic Dirac cone. The absolute value of the component 
of the eigenvector is independent of 0, i.e. |d Q7 (0)| = l/\/2. 
The eigenvector of the latter case may correspond to the 
component of A and A' of the former case, in which the 
0-dependence is small as shown in Figure [3] Thus it is 
Jbund that the existence of the B and C components of 
the present case is the characteristic of the Dirac particle 
of organic conductor, a-(BEDT-TTF) 2 I 3 salt, in which 
there are four sites in the unit cell. Further the vanishing 
of the B and C components at a certain value of seems 
to be general and exist for other case of transfer energies 
as discussed in f36l 



4 The effect of the anion potential 

Next, we examine the anion potentials, Ia(= h = I2), 
Ib{= I3), an d Ic(= I4) perturbatively. Although the charge 
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Fig. 5. Ellipsoidal curve obtained from the condition, £ = u), 
for the Dirac cone. 6 is the angle between k — feo and fe^-axis, 
where feo is the point on the ellipsoidal curve. The inset denotes 
a tilted Dirac cone where uj = corresponds to the apex, i.e., 
the Dirac point. 

disproportionation depends on the anion potential [38], 
quantities of (n a ) are fixed to (ni) = (n 2 ) = 1.460, (n 3 ) = 
1.369, and (714) — 1.711, which is obtained by the condi- 
tion, 1 = 0. Such a treatment can be justified for the qual- 
itative discussion of the location of feo since the symmetry 
relation, (ni) = (n 2 ), remains due to I\ = I 2 . 

When the perturbative Hamiltonian is given as 



V(k) = \k ,a) 
eq. (fTTj) is rewritten as 



df 
k -dk {ko) 



{k o ,0\, (25) 



[-K<C" + k ■ ««" + Jc<"W"(> (fc) = AE c (k)a cc (k) 

C" = l,2 

(26) 



with 



£ 

77=3,4 



(fc ■ V( v )(k ■ Vnc) 
L ( fe o) 



"cc = vtf - 



j w = 1 



where the matrix elements, 7 77 ' (7 = C? 7 ? = 1) 
defined by 



(27) 
,4) is 

(28) 



Discarding the first term of eq. (|2"6"|) for simplicity, i.e. 
K 11 ' — > 0, the energy bands are calculated as 



AE 1)2 (fc) =k-u + J ± 



■ «, + J P ) 2 , (29) 



Table 1. Quantities of coefficients, K" x , K" y (a = A,B,C, 
P= 0,1,2,3). 



J px 


p = 


p = l 


p = 2 


p = 3 


A 


-0.03 


0.00 


0.01 


0.12 


B 


0.08 


0.00 


0.14 


-0.07 


C 


-0.04 


0.00 


-0.15 


-0.04 


K a 
py 


p = 


p = l 


p = 2 


p = 3 


A 


0.01 


0.00 


0.07 


0.00 


B 


0.01 


0.00 


0.02 


0.00 


C 


0.00 


0.00 


-0.10 


0.00 



where u = («n+i*22)/2, «i = Re(Mi 2 ), u 2 = -Im(ui 2 ), 
u 3 = (un - u 22 )/2, J = (Jn + J 22 )/2, Ji = Re(Ji 2 ), 
J 2 = — Im(Ji 2 ), and J3 = (Jn — J 22 )/2. Compared with 
the previous section, the velocity is changed by the anion 
potential. Actually, for the bases chosen as eqs. (fTO]) and 
(HHJ), components of u p (p = 0, • • • ,3) are given by 



l pX 



J px 



K±I A 



■K?J B 



K c T 

px^si. ' ' px C 5 



py^B py^C 



(30) 



where coefficients of 



iT^ are given in Table CD The 



velocity vector as a function of Ic is drawn in Figure [H 
where —0.2 < Ic < 0.2 and the other anion potentials are 
set as Ia = Ib = 0. In the presence of J p , the position of 
the contact point changes from that of the state with I = 
0. The quantities J p (p = 0, 1, 2, 3), which are expanded up 
to a quadratic term of the anion potential, are described 

by 



- ab 



tCj ,tAAt2 



T BB t1 
L p 1 B 



T CC t2 
L p 1 C 



L a I a Ib 



L p c IbIc 



L^ A IcIa , 



(31) 



where L" is given in Table (2] The Tp-dependences of J\ , 
J 2 , J3 (I a = Ib = 0) are shown in the inset of Figured] 
Note that Ji and J 2 are negligibly small compared with J3 
(and also for case of the /^-dependence with Ia = Ic = 
0). It should be noted that the effect of Jo being the same 
order of J3 does not contribute to the shape of the Dirac 
cone. 

Using the above results, the energy difference between 
AEi(fe) and AE 2 (k) is approximately given by 

AEa(fc) - AE 2 {k) 

~ (u 2x k x + u 2y k y ) 2 + (u 3x k x + u 3y k y + J3) 2 . (32) 

Then, the contact point, which is given by AE\ — AE 2 = 0, 
is obtained at 



J3U 2y 



^0x 



U2xU 3 y - U 3x U 2 y 



{ 0y 



-J3U2X 



U 2x U 3 y - U 3x U 2 y 



(33) 



The trajectories of the contact point calculated from both 
eqs. ([6]) (circle) and (|33p (line) are drawn in the main part 
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Table 2. Quantities of coefficients, L™ (a 
A, B, C, AA, BB, CC, AB, BC, CA). 





p — U 


P - 1 


p — 2, 


p — 6 


L p 


0.47 


0.00 


-O.Of 


0.07 


I B 

L p 


0.30 


0.00 


0.01 


-0.30 


tF 
l p 


0.23 


0.00 


0.00 


0.23 


r AA 
L P 


0.86 


0.00 


-0.26 


0.01 


r BB 
L P 


0.42 


0.00 


0.01 


-0.42 


rCC 
L P 


0.44 


0.00 


0.00 


0.44 


t AB 
L P 


-0.85 


0.00 


0.25 


0.85 


r BC 
U P 


0.00 


0.00 


-0.27 


0.00 


r CA 


-0.87 


0.00 


0.27 


-0.87 



of Figure [51 In case of the exact calculation, it moves from 
Y point to M point as shown in Ref. [58] . The trajectory 
obtained from eq. ([55)) agrees qualitatively with that of 
eq. for the small potential. Thus, one finds the ZGS is 
robust against the anion potential. 



r x 




Fig. 6. Trajectory of the contact point in the quarter part of 
the Brillouin zone (0 < k x < n, —n < k y < 0). Filled (open) 
circle denotes the position of feo calculated from eq. ([6} by the 
condition I A = 0, -0.3 < I B < 0.3, I c = (I A = 0, I B = 
0,-0.3 < Ic < 0.3). The corresponding result obtained from 
eq. (|33|) is drawn by solid (dashed) line. The inset shows the 
quantities Ji, J2, and J3 as functions of Ic- 



5 Local density of states and magnetic 
properties 

By using the energy band and eigenvector (eq. ©), the 
density of states of the a site in a unit cell is calculated 
as 

= jf E K-r( k )\ 2 ^ ~ er(*))> (34) 

where a = 1(A), 2(A'), 3(B), and 4(C) denote the respec- 
tive sites. The total value of the density of states is given 
by J deJ2a=i D a (e) — 4 and e = denotes the chemical 
potential. Figure [7] depicts D a (e) close to the Fermi en- 
ergy. The local density of states, D a (e), for small e is the 
largest at C site corresponding to the electron rich site, 
and is the smallest at B site corresponding to the hole (or 
charge) rich site [3D]. We notice that the above result is 
consistent with that of the experiment [33] . 

The relative magnitude of D a (e) is understood as fol- 
lows. From Sj3j we obtain that d aj (k) depends on only 
9(= tan _1 [fcy/fc :E ]) in the vicinity of the Dirac point, and 
that the energy band is described by eq. ((24|) . Then D a (e) 
(e > 0) can be rewritten as 

I I />27T />OC 

D a {e) = lr\ d$ drr\d al (0)\ 2 S(e-r/F(9)) 
^ Jo Jo 

= K &9\d al (9)\\F{e)f, (35) 

and 

I>( £ ) = fi M{F(0))\ (36) 




Fig. 7. Local density of states D a (e) per unit cell for A and A' 
sites (dotted line), B site (dashed line), and C site (dot-dashed 
line) , which are calculated from eq. (|34|l . Solid line is the total 
density of states. The parameters of interactions are taken as 
U = 0.4, V a = 0.17 and V b = 0.05. 



in the low energy region. The integration of eq. (I35p is 
dominant for 9^0, since F(0) is larger than F(tt) as 
shown in Figure [5] Further, from Figure [3] (a), we see 
that the component of the eigenvector satisfies a condi- 
tion |du(0)| > \d n (0)\ = \d 21 (0)\ > |d 3 i(0)|. These facts 
give a conclusion that the contribution of the C site (elec- 
tron rich site) to the density of state is larger than that 
of the B site (hole rich site). It is also found, within the 
present scheme of eq. (fTB]) . that the density of states of 
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e < is the same as that of e > 0, from the symmetry 
of the energy band, AEi(k) = -AE 2 (-k) (see eq. (fl5]>). 
The underlying symmetry is due to the transformation 
property of the effective Hamiltonian under space inver- 
sion [17]. We note that eq. ([M]) is the same as \e\/ (2ir 2 v F 2 ) 
obtained in Ref. [17] . The qunatity v* F is the renormalized 
velocity, and is estimated as v* F = 0.0097 for the choice 
of the present parameters where the corresponding four 
effective velocities in Ref. [17] are given as wo x = —0.039, 
w 0y = 0.008, w x = 0.052, and w y = 0.044. For small \e\, 
the total density of states shown in Figure [7] coincides well 
with both cq. ([36]) and \e\/(2n 2 v* F 2 ). 

Finally, the magnetic properties are examined. The 
susceptibility is calculated on the basis of four sites in a 
unit cell. The susceptibility corresponding to the response 
between the a site and the (3 site is described as, 

X%{Q^) = 

T7 d aj (k + q)d*p 7 (k + q)dp Y (k)d* aY (k) 



N 

£ 7 (fe + q) — £ 7 / (fe) — lj — id ' 



(37) 



where i<5 (6 > 0) is a infinitesimally small imaginary part. 
Using equation ([3"7]) . the magnetic responses of both the 
local susceptibility and the local NMR relaxation rate [41] 
14121 at the a site are obtained as 



Xa = ^Re( X ^(O,0)) = f deD a {e) (- 

J — OG \ 



ds 



— =T lim 



ttT 



X 



de(D a (s)Y 



df 

ds 



(38) 



(39) 



respectively. 

The local magnetic susceptibility and the local NMR 
relaxation rate for each site are shown as a function of 
temperature in Figures [8] and [9] respectively. The suscep- 
tibility shows a relation \c > XA — XA' > Xb, which 
agrees qualitatively with that of the experiment [34] . The 
open circles denotes the total one, which is calculated us- 
ing the temperature dependence of the chemical potential 
and (n a ) (eqs. ©-([8])) and the transfer energy given by 



t L (T) 



t £ (300K)-t L (8K) 
292K 



(T-300K) + < L (300K). 



(40) 



[eV _1 ] 6 







o 


/a_i- a '_i_R_i_r 

c 




"""" A, A' 




B 



0.005 



0.01 



T 
[eV] 



0.015 



Fig. 8. Temperature dependence of \ a (a — A,A',B,C). 
Dotted line, dashed line, and dot-dashed line correspond to 
A (A'), B, and C sites, respectively. The total value of \a is 
plotted by the solid line where the open circle is calculated 
using eq. (J40j . 



([39]) at low temperature are calculated as Xa = 2(\n2)K a T 
and (l/Ti) Q = (ir 3 /3)K*T 3 where D(e) = K a \s\ from eq. 
([3"5]. where K A ~ 137.1, K B ~ 38.9 and K c ^ 212.9. 



The data of t L (300K) and t L (8K) (L = al, ■ ■ ■ , a3, 61, • • • 
al',a3',a4') are given in Ref. [3H]- The deviation is visi- 
ble for T > 0.005 eV. In Figures E and M reflecting the 
density of states, Xi an d (1/^1)4 (C site) are the largest 
among those of the four sites. It is found that, at low tem- 
perature, Xa and (1/Ti) a are proportional to T and T 3 , 
respectively, due to D a {e) cx |e|. Actually, eqs. ([38]) and 



6 Summary and discussion 

We examined the ZGS close to the Dirac point, and local 
magnetic properties at low temperature. Each electron in 
the unit cell has respective role to form a Dirac particle 
which is different from that of the graphene. The C and 
B sites give the largest and smallest density of states re- 
spectively, and are peculiar for the present salt a-(BEDT- 
TTF)2l3 while the A and A' sites have a common feature 
with the graphene. The temperature dependence of mag- 
netic susceptibility and NMR relaxation rate exhibit the 
power law at low temperature, which can be interpreted 
in terms of the density of states close to the Fermi surface. 

Here we note that the detail of the ZGS depend on the 
choice of parameters as transfer energy, interaction and 
pressure. There are two kinds of conditions for the exis- 
tence of the ZGS. First one is the contact point, which 
is satisfied for most of a- type BEDT-TTF salt, and may 
arise from a property of an inversion symmetry of the four 
sites in the unit cell. The second one is that the contact 
point exists on the Fermi energy (i.e, the disappearance of 
hole or electron pockets). The ZGS has been obtained in 
the following cases, (i) Using the data of the X-ray exper- 
iment at room temperature [6] and U = 0.4, V a = 0.17, 
V a — 0.05, the ZGS occurs for P a > 4.3 kbar, whereas the 
64fharge ordered insulating (metallic) state is obtained for 
kbar < P a < 3.3 kbar (3.3 kbar< P a < 4.3 kbar) (3|. 
(ii) In the absence of interaction, the ZGS (metallic) state 
is obtained for P a > 3 kbar (3 kbar< P a ) [5]. (hi) For 
the small magnitude of interaction, the ZGS is obtained 
due to reducing the pockets while the large interaction de- 
stroys the ZGS due to the formation of the charge ordered 
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(l/7i) 
[eV _1 ] 




D«(e) 



Fig. 9. Temperature dependence of (1/Ti) a - Dotted line, 
dashed line, and dot-dashed line represent (1/Xi) a for A (A'), 
B, and C sites, respectively. The total value of Xq is plotted 
by the solid line where the open circle is calculated by using 
eq. @0]|. 



(CO) state, (iv) Further, the ZGS appears under higher 
hydrostatic pressure (12.5 kbar < P) [43] , 

We note the following relation between the charge dis- 
proportionation and the magnetic susceptibility in the ZGS. 
Using parameters at P a —10 kbar and room tempera- 
ture [6], the charge disproportionation is given by (nc) > 
(ua) = (tia') > (n B ) while the susceptibility (oc T) is ob- 
tained as xa(— XA') > Xc > Xb- This is contrast to the 
case for transfer energies calculated from the first princi- 
ple calculation at 8K and at ambient pressure [3D, where 
(nc) > (tia) = (nA>) > (n B ) and X c > Xa = Xa> > Xb- 
For the present choice of the transfer energy [36] , the ZGS 
(CO) state is obtained for x < 1.03 (x > 1.03) with 
U = 0.4x, V a — 0.17a; and V& = 0.05a;. Thus, we propose 
the following variation of interaction in order to explain 
the experiment that the CO state is obtained at ambi- 
ent pressure and the ZGS is obtained for P a > 4.3 kbar. 
The parameter of interaction at P a = kbar is taken as 
x > 1.03 and that of P a > 4.3 kbar is expected as x < 1.03, 
since the effect of pressure increases the band width and 
decrease the relative magnitude of interactions. 

Finally we discuss the state, which is derived from the 
transfer energies obtained from the experimental data at 
T = 8 K and at ambient pressure [6]. In this case, the 
ZGS emerges even at ambient pressure. Figure [TO] shows 
the density of states under the ambient pressure where 
the parameters of interactions are the same as those of 
Figure 2. The charge disproportionation is obtained as 
(nc) > (tia) — (riA') > (ib) which is qualitatively the 
same as the first principle calculation. [35] However, com- 
pared with Figure!?] the difference of each of D a (e) is very 
small in the vicinity of the Fermi energy. This implies that 
the similarity between Xa (corresponding to D a (e) ) and 
(n a ) is not necessary for the existence of the ZGS. We also 
note that the node found in Figure[3jalso exists in the case 
of Figure [TO1 although the location of the node is different. 
Thus the role of respective sites seems to be a common 




0.02 



Fig. 10. Density of states at the ambient pressure and at 
T = 8K where interactions are the same as Figure [7] 



feature for the ZGS of the system with four sites in the 
unit cell. Further we comment on a large Van Hove sin- 
gularities, which exist for e > in Figure [71 and for e < 
in Figure [TO1 These singularities give opposite results for 
temperature effect, e.g. with increasing temperature, the 
chemical potential decreases for Figure [7] but increases for 
Figure [TO] with increasing temperature. Thus the case of 
Figure [TO] needs further consideration in order to explain 
a fact [TO] that the behavior of the Hall coefficient sug- 
gests the monotonic decrease of chemical potential with 
increasing temperature. 
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Appendix: A constant energy curve for the 
general tilted Weyl equation 

In this appendix, we examine the relation between the 
linear dispersion calculated from the tilted Weyl equa- 
tion and the ellipse as shown in Figure [5] When we put 
A+(fc) = u> in eq. ([TB"]) . we obtain the general equation of 
the ellipse centered on (Xq,Yq) as follows: 



V 2 



"0.T 



(k 3 



R 2 

V xy ~ «0^0j/ 

R 2 



X n ) 2 



V: 



'0y 



R 2 



Yn) 2 



(k X -X Q )(k y -Y ) = 1, 



(41) 
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where 

R = LJ 



y2y2 _ yi 
y x y y v xy 



-,1/2 



(V 2 - wgJOtf - v 2 y ) - (V x % - v 0x v 0y y 



(42) 



The coordinates of the center (Xq,Yq) are given in the 
form: 



= 

w (V 2 

^0 



VOyV* 



y v °* V v 



«ta)(ty " V 0y) ~ ( V xy ~ V 0x V 0y ) 2 



voxV^y - VQyV 2 



» (Yi - vl x ){V* - v$ y ) - (V 2 y - v 0x v 0y y 



(43) 



Carrying out the coordinate rotation, 



cos <p sin < 
— sin 6 cos ( 



4>\ ( k x - Xq 



ky - Y 



• tan 



7T 

2 



2{V xy ~ VQxVQy) 



(V 2 



(K 2 



- "0y> J 

vL < V 2 



tan 



"Vx " y 
2 (Vxy ~ VQxVQy) 



v Qy) 



( V x 



' 0.T 



> VI 



v Qy) 



(44) 



we obtain the standard form of the ellipse (Fig. [5]) given 

by 



a 2 b 2 ' 
where the long and short axes are given as 

V2R 



(45) 



V 2 



1/2 ' 



V2R 



V 2 



1/2 



D 



[(V 2 - V y 2 ) - (v 2 0x - v 2 0y )} 2 + A[V 2 y - v 0x v 0y } 2 . (46) 
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